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ABSTRACT 

The low-mass end of the stellar Initial Mass Function (IMF) is constrained by focusing on 
the baryon-dominated central regions of strong lensing galaxies. We study in this letter the 
Einstein Cross (Q2237+0305), a z=0.04 barred galaxy whose bulge acts as lens on a back- 
ground quasar. The positions of the four quasar images constrain the surface mass density on 
the lens plane, whereas the surface brightness (H-band NICMOS/HST imaging) along with 
deep spectroscopy of the lens (VLT/FORS1) allow us to constrain the stellar mass content, for 
a range of IMFs. We find that a classical single power law (Salpeter IMF) predicts more stel- 
lar mass than the observed lensing estimates. This result is confirmed at the 99% confidence 
level, and is robust to systematic effects due to the choice of population synthesis models, the 
presence of dust, or the complex disk/bulge population mix. Our non-parametric methodology 
is more robust than kinematic estimates, as we do not need to make any assumptions about 
the dynamical state of the galaxy or its decomposition into bulge and disk. Over a range of 
low-mass power law slopes (with Salpeter being T = + 1 .35) we find that at a 90% confidence 
level, slopes with T > are ruled out. 

Key words: galaxies: individual: 2237+0305 - galaxies: stellar content - stars: luminosity 
function, mass function - gravitational lensing. 



1 INTRODUCTION content. Furthermore, changes in the low-mass end of the IMF do 

not translate into significant changes of the photo-spectroscopic ob- 

The stellar Initial Mass Function (IMF) is defined as the start- . . . . , „ , ,. ....... , 

. . servables, mainly because of the small contribution from low-mass 

ing mass distribution of a stellar population. A large number of . , . , , . - . TI , 

° r 1 ° stars to the net luminosity budget ot a galaxy. However, changes in 

studies have targeted its shape and possible universality (see e.g. , , - „„ . , , , 

I ii 1 i " — ii n f —ft -3 rr — 7— rii . i the low-mass end ot the IMF result in measurable changes of the 

Kroupa 2002; Larson 2006; Kroupa 2007; Bastian et al. 2010). It „ „ . . .. , . . . _ , * . . 
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has important implications in the formation and evolution of galax- 



total stellar mass of a given population. Therefore, any attempt at 

constraining the low-mass end of the IMF requires an independent 

ies as it affects the chemical enrichment process and determines r . r , rz; rr - : — ~r\ 

r . I — 1 measure ot the stellar mass content ot the galaxy. Cappellan et al. 

the stellar mass content of galaxies. The seminal work of Salpeter ] n „„i , , ,. , „ , ro«TTr.^T 

i — — i (2006) compared the kinematic data from a number of SAURON 

( 1955) assumed a single power law distribution between the upper , . , . ... , . , , , , , , , . 

- 1 ' early-type galaxies with dynamical models, and concluded that a 



(~1OO-12OM ) and lower mass cutoffs (0.07M©). More recently, 

the "local IMF" has been constrained in more detail, with a turnover ........ . . 

ing gravitational lensi ng in the inner, bary on-dominated, regions 



single power law for the IMF (a la Salpeter) was ruled out. Us- 



i -ii -i ing gravitational lensing in tne inner, baryon-aominated, regions 

at lower masses (Andersen et al. 2008), and a characteristic mass r . , , . rrr — TTIAorxiN i j .i t o i . 

i 1 ii I ot early-type galaxies, Ferreras et al. ( 2008) showed that a Salpeter 

aroun ®" arson ll W- IMF produced too much stellar mass compared to the lensing mass 

Constraining the IMF outside our Galaxy is a challenging .... T . , .. .. , 

° J 00 within the same region. In this paper we extend the gravitational 

task. Data from unresolved stellar populations are very degener- . . , . ... , „. , 

r r lensing analysis and the modelling work, to give a more quan- 

ate with respect to their properties such as age, metallicity or dust ... .. . . . . , . „„ „. . , . 

r r r ° J titative constraint on the low-mass end of the IMF. We study in 

this letter a nearb y strong lens, the Einstein Cross (Q2237+0305 



tms letter a nearby strong lens, tne tinstein Cross /+UoUj 
iHuchra et alJI 19851) . combining high resolution NIR imaging with 
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Figure 1. Cumulative lensing mass profile (shaded region) and luminosity 
(dots). 90% confidence levels shown. Notice the butterfly-shape of the un- 
certainty for lensing mass, which is minimal at the position on the lens plane 
of the images of the background source. The inset shows a contour map of 
the F160W HST/NICMOS image. At the lens position 1 arcsec maps into 
0.77 kpc. The arrow points North, the small segment from the base of the 
arrow points East. 

HST/NICMOS and a deep spectrum of the lensing galaxy taken 
with the VLT, using the FORS 1 instrument, to constrain the prop- 
erties of the stellar populations. 

The reason for choosing this lens is the exceptionally low red- 
shift (z=0.04) of the lensing galaxy. This has two advantages. First, 
it helps increase the signal-to-noise. Second, and even more impor- 
tant, the low lens redshift makes a strong lens from the innermost 
and densest part of the galaxy. The mass-to-light ratio is close to 
unity in the lensing region (see Figures 1 and 2 below). In other 
words, the lens is dominated by stars, and hence particularly well- 
suited to constraining the IMF. In contrast, typical galaxy lenses 
contain a significant dark-matter contribution. 

A di fferent app r oach u sing lensing galax i es has been recently 
taken by iTreu et all d2010h and lAuger et all d20ld) . Rather than 
choosing a single lens with minimal dark matter and maximal 
signal-to-noise, they choose a large sample of lenses, and attempt to 
subtract off the dark matter by fitting dark-matter profiles based on 
simulations. Their results fa vour a Salpeter I MF. However, a caveat 
in the method (as noted bv lTreu et al]|201Cl) is that if the assumed 
dark-matter profile is slightly too shallow, the stellar mass will be 
systematically overestimated, and a "heavy" IMF derived. If this is 
in fact the case, that wo uld explain the apparent contradiction with 
the lFerreras et al.l d2008h result, which excludes a Salpeter IMF for 
a subset of the same galaxies. 



2 CONSTRAINING THE IMF 

In order to constrain the low-mass end of the IMF, we need to com- 
bine stellar population synthesis with an independent estimate of 
total mass in a system whose total mass content is preferrably (but 
not necessarily) dominated by the stellar mass. We note that even 
though baryonic, non-stellar matter (e.g. gas and dust) will con- 
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Figure 2. Predicted H-band mass to light ratio combining the analysis of 
the lensing data and the luminosity profile of the lens. The dot corresponds 
to the value chosen in this paper to constrain the stellar M/L, with mean and 
standard deviation as labelled. 



tribute to the mass budget, the constraints on the IMF imposed in 
this letter go in the sense that rejected choices of the IMF predict 
more stellar mass than the total mass determined by the lens geom- 
etry. We use strong gravitational lensing as the independent esti- 
mator of mass. Kinematics-based methods have been used recently 
to quantify the total mass in the inner regio ns of galaxies (see e.g. 
ICappellari et alj|2006t ICoccato et ai]|2009h . However, such meth- 
ods rely on a particular definition of the dynamical system. 

The lens is mo d elled using the PixeLens method 
dSaha & Williams! |2004 IColesI \200A whereby the projected 
mass distribution is reconstructed as a pixell ated mass map, which 
is not limited to any parametric form (cf . iTrott et all l2010h but 
can be any non-negative distribution subject to the following prior 
conditions, (i) The mass map is symmetric under a 180° rotation 
about the centre of brightness, (ii) The local density gradient must 
point no more than 45° away from the centre, which is assumed 
to coincide with the light peak, (iii) No pixel, except the central 
one, is allowed to be more than twice the mean of its neighbours, 
(iv) The circular average of the projected density, falls off as 
i?~ 1//2 or steeper. In other words, the effective three-dimensional 
mass profile should be oc r~ 3 ' 2 or steeper. 

The above requirements do not give a unique mass map, 
because solutions of the len s equation are highly n on-unique 
dFalco et alj|l985l ; ISahal kood : lliesenborgs et alj|2008h . Accord- 
ingly, the possible "model space" is sampled by an ensemble of 
300 models. A useful property of the ensemble, due to the linear 
nature of the observational and prior constraints, is that the ensem- 
ble average is also a valid model; for comparisons of ensemble av- 
erages with indiv i dual m odels in PixeLens, see Figures 7 and 8 in 
ISaha& Williams] ( 120061) . 

F or this study, we select the Einstein Cross jHuchra et al.l 
1985), where the bulge of a barred Sab galaxy at z=0.0394 acts as a 
lens, generating a quadruple image of a background z= 1.695 QSO. 
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Figure 3. The 68, 95 and 99% confidence levels in metallicity and average 
age are shown for the r models described in the text, using the Pegase-HR 
(High Resolution) library. The greyscale corresponds to an IMF with slope 
r = +1.35 (i.e. Salpeter), whereas the lines represent the other extreme 
explored in this paper (r = —1.4). 

Our reference image is taken from the CASTLES survejQ. It is a 
384 second exposure with NIC2/HST through passband F160W. 
The lensing analysis only uses the redshifts of the QSO and of the 
lens, along with the QSO image positions to constrain the surface 
mass density on the lens plane. The NICMOS image is also used to 
recon struct the surfac e brightness profile of the lens. We use GAL- 
FIT jPeng et al.l2 002) in or der to remove the unresolved images of 
the background QSO (see iLeier et alJkoid for details). Figure [TJ 
shows the cumulative mass and H-band luminosity profiles of the 
lens. The uncertainties in lensing mass and luminosity are given by 
the shaded region and errorbars, respectively. Notice the lens ge- 
ometry allows us to constrain in detail the cumulative mass within 
R ^ 0.75 kpc, roughly around one effective radius (top axis). 

The range of lens models has a characteristic butterfly shape, 
which can be understood as follows. For a pure Einstein ring, the 
enclosed mass M(< R) is perfectly constrained at the Einstein 
ring, but degenerate at any other radius. Now, the prior allows mass 
profiles ranging from R~ x ^ 2 to a point mass, hence M(< R) can 
range from R 3 ^ 2 to flat. The maximal and minimal values of steep- 
ness, combined with the Einstein-radius constraint, shape the al- 
lowed range of models. In practice, because real lenses are not 
pure Einstein rings, the range of models is not quite so extreme, 
but retains the butterfly shape. The inset shows the H-band frame 
of the lens and quasar images along with their standard identifica- 
tions. The mass and luminosity profiles are compared to obtain the 
H-band mass-to-light ratio (Th), as shown in figure [2] including 
the uncertainty as a shaded region. The point where the errorbar is 
smallest is chosen to constrain the low-mass end of the IMF. The 
figure shows that within this region, the lensing galaxy has an av- 
erage T H = 1.19 ± 0.16 (90% confidence level). 

Our next step involves the construction of a large number of 
composite stellar populations, following an exponentially decay- 
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Figure 4. Same as figure [3] comparing two different sets of models, both 
assuming a Salpeter IMF. The greyscale corresponds to Pegase-HR, and the 
lines are the predictions of the BC07 models (see text for details). 



ing star formation history. Hence, each model is characterised by 
the formation epoch, decay timescale and metallicity. We also need 
to consider different sets of models for a number of choices of the 
IMF. Notice these are composite models and not simple stellar pop- 
ulations (SSPs), therefore giving more robust mass-weighted val- 
ues. SSP-determined values are easily affected by the presence of 
recent star formation, even with small amounts of young stars, re- 
sulting in - some times significantly - lower M/L ratios (see e.g. 
iRogers et alj20ld) . 

In previous work jFerreraset J120 05. 2008; Leier et alj20ich 
we used broadband photometry to constrain the mass-to-light ra- 
tio of the stellar populations. We showed that acceptable uncer- 
tainties can be obtained on M/L even though one should not use 
broadband photometry alone to constrain the age and metallic - 
ity of the unresolved populations (see e.g. iGallazzi & Belill2009l) . 
However, in order to minimise the uncertainties in the estimates 
of M/L, we use deep spectroscopic observations of the Einstein 
Cross from VLT/FORS1 wi th grism G300V, givin g a spectral reso- 
lution of R~ 400 at 5900A (Eige nbrod etalj2008l) . The spectra are 
taken at two slit positions: mask 1 (tEXP=H-70h) at PA= —56.5°, 
passing through images A and D; and mask 2 (tEXP=lL25h) at 
PA= —68.5°, passing through images B and C (see figure[T). The 
2D spectra are separated into a point source - the quasar spectrum 
- and an extended source - the lens. The method leaves no trace 
of the emission line spectrum of the quasar in the len s spectrum, 
which can be seen in figure 7 of Eigenbrod et al. ( 200^). 

We measure the equivalent widths (EWs) of a number of age 
and metallicity sensitive spectral features: Balmer H/3, H7, H5, 
Mgb, Fe5270, Fe5335, G4300 and the 4000A break. The EWs are 
determ ined using a recent definition of the continuum (Rog ers et al.1 
l20ld) that reduces the contamination in EW estimates from neigh- 
bouring lines, improving on the discrimination between age and 
metallicity. 

We also have at our disposal the HST/NICMOS (F160W) and 
HST/WFPC2 (F555W and F814W) images. Note we already use 
the F160W image to determine the maximum Th from stars. How- 
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ever, we can combine these three bands to generate two colours. At 
the effective radius, the colours are V — I = 1.28 and / — H = 
1.89 (Using Vega as zero-colour reference). The la photometric 
uncertainties stay around 0.05 mag in both colours. The colour 
gradient is very shallow: over the R< 2R e range, the gradients 
AColour/A(i?/i? e ) are -0.003 and 0.05 in V - 1 and / - H, re- 
spectively. The flat colour profile confirms the homogeneity of the 
stellar populations in this lens. The colours help further constrain 
the stellar populations, although we emphasize that the inclusion of 
the colours does not improve significantly the estimates based on 
EWs, as expected. 

From the available stellar population synthesis models, 
we chose the high spectral resolution version of Pegase 
dLeBorgneetaT]l2004h . These models allow the user to recalcu- 
late the spectra of simple stellar populations for an arbitrary choice 
of IMF. Our goal is to constrain the shape of the low-mass end of 
the IMF. Hence, we decided to model the IMF as two power laws, 
defined as follows: 



dK 



dlogM 



M" 



M" 



0.1 < M/M < 1 
1 < M/M < 100 



(1) 



The above definition fo llows Kroupal |2007|) . although we select 
the Salpeter IMF slope (Salp etej|l95a) at the high mass end and 
choose IMq for the position of the knee. For masses below this 
knee, we parameterise the IMF as a power law with index T. We 
explore a wide range of values, from T = —1.4 to +1.4 (where 
a Salpeter IMF corresponds to T = +1.35). For each choice of 
r, we generate a set of base models (simple populations) that are 
combined into composite r-models as described above, covering 
a grid of 48 x 48 x 48 models over a range of formation epoch 
(2 zfor ^ 10); decay timescale (—1 log(r/lGyr) +1) 
and metallicity (—1 ^ [m/H] ^ +0.3). Hence, our methodology 
takes into account a spread in the age distribution, an important is- 
sue when dealing w ith the stellar populations of galaxies (see e.g. 
Roger s~et alj|201oh . Each choice of star formation history gives a 
set of EWs and colours that are compared with the observations 
via a standard \ 2 -defined likelihood to infer the probability distri- 
bution function. The synthetic spectra are broadened to the spectral 
resolution of the observations, taking also into acc ount the observed 
veloc ity dispersion of the lens, a = 166 + 2 km/s dvan de Ven et alj 
2008). Also notice that the constraint on the M/L ratio shown in fig- 
ure|2]corresponds to the observed frame H-band. The models take 
into account the minor K-correction needed, given the redshift of 
the lens. 



3 THE SYSTEMATICS OF STELLAR MAS ESTIMATES 

It is important to understand the systematic effects caused by the 
choice of population synthesis models in the estimates of stellar 
mass. Firstly, we consider the effect that a different choice of slope 
would have on the predicted average age and metallicity of the 
lens galaxy. Figure [3] shows the probability distribution of the aver- 
age age and metallicity between the two extreme choices of IMF: 
namely a Salpeter IMF (T — +1.35, greyscale) and a strong sup- 
pression of low-mass stars (r = —1.4, lines). One can see that 
within error bars, the average age and metallicity are unchanged, 
corresponding to an old ( ^ 8 Gyr) population. 

Different population synthesis models rely on their choices of 
stellar evolution prescriptions, isochrones, and stellar libraries. We 
quantify this effect in the following two figures. Figure [4] shows 
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Figure 5. Comparison of the H-band stellar mass-to-light ratio predicted 
by a number of population synthesis models. For ease of comparison, we 
consider here only simple stellar populations with solar metallicit y and a 
Salpeter IMF. The models are Pegase (black solid line, Le Borg ne et alj 
|2004|): two different flavours of the GAL AXEV models: BC03 (blue dotted 
line lBruzual & Charloj2 003 ) ; and the recent, unpublished 2009 models in- 
cluding the MILES spectral library (orange shor t dashed lines, private com- 
munication, and Sanchez-Blazquez et al. 2006); Vazdekis/MILES models 
v.9 (l ong da shed lines lVazd^kiset*all20*lol) "and lMarastonl (grey dot-dashed 
lines. 120051) . The bottom panel shows the ratio with respect to the fiducial 
Pegase-HR models used in this paper. 

the age and metallicity probability distribution for a Salpeter IMF 
using two independent models: Pegase (grayscale) and Bruzual & 
Chariot (2003; lines). The latter tends to give slightly older ages, 
although the confidence levels overlap, eliminating the chance of a 
systematic offset caused by a different age or metallicity. 

To emphasize this point, figure [5] illustrates the predictions 
in Th for a range of Simple Stellar Populations of solar metal- 
licity - as suggested by the best fits shown in the previous fig- 
ures. Various population synthesis models are shown, as labelled. 
Hence, for the age range considered to be a best fit ( £ 8 Gyr), the 
systematic difference between these models cannot be larger than 
~ 10%. Another possible systematic source wo uld be the pres- 
ence o f dust, mainly in the disc of the lens galaxy. lEigenbrod et al.1 
(2008) showed that the U-band differential extinction caused by 
the lens on the background quasar is in the range 0.1 -0.3 mag. Us- 
ing a standard extinction law (e.g. iFitzpatrick 1 19991) , we estimate 
an attenuation no larger than about 0.05 mag in the H band. Hence, 
the net systematic uncertainty on the determination of the stellar 
masses is ~ 15 — 20%. 



4 RESULTS 

Figure [6] shows the constraint on the stellar Th from the photo- 
spectroscopic data, as a function of IMF low-mass slope. For ref- 
erence, the upper limit to the Th from lensing is shown as a hor- 
izontal line along with the 90% confidence level (dashed lines). A 
Salpeter IMF (r = +1.35) clearly gives a stellar Th that is too 
high with respect to the lensing estimate. The best fits are shown 
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Figure 6. Top: Constraints at the 90% confidence level from the observed 
equivalent widths of Balmer lines (H/3, H7, H5); the metal sensitive [MgFe] 
index; and the V — I and I — H colours. The result is shown as a func- 
tion of the low-mass slope of the IMF. The shaded region corresponds to 
the Tg given by the Pegase-HR models used in this paper. The horizontal 
lines give the 90% confidence level for the predicted M/L using lensing data 
and the observed photometric profile of the lens (see figure [5). The solid 
(long dashed) curved lines correspond to slit 1 and 2, respectively, with the 
shaded region extending over the 90% confidence level for both slits 1 and 
2. Bottom: Cumulative joint probability of the IMF slope, given the lensing 
estimate of T h and the population synthesis constraints. A Salpeter IMF 
(r = +1.35 is ruled out at more than the 99% confidence level. At a 90% 
level, slopes T > are rejected as well. 

as curved solid and long-dashed lines for the analysis of spectra 
from masks 1 and 2, respectively. The result is compatible within 
the 90% confidence level (grey shaded area). Furthermore, our data 
can constrain in more detail the low-mass end of the IMF. The 
lower panel of figure [6] shows the cumulative probability distribu- 
tion function, P(> T). At a 90% confidence level, slopes Y > are 
ruled out. 

Our modelling has taken an "integrated approach" to the pres- 
ence of disk and bulge. Indeed, both the spectroscopic and the pho- 
tometric data have not been separated into a bulge and a disk com- 
ponent. Hence, the modelling of the data corresponds to an "av- 
erage" stellar population which would be dominated by the old 
component, given that most of the light originates from the central 
region of the galaxy, a point confirmed by the old ages and high 
metallicities predicted by the modelling (see figures[3]and[4]l. The 
photometric analysis is done in a non-parametric way, hence we 
compute stellar masses from this "composite population" of disk 
and bulge stars in a consistent way. Were we to remove the pres- 
ence of the disk population, the older bulge stars would result into 
higher Th therefore making a Salpeter IMF even more incompat- 
ible with the observations. Finally, a significa nt presence of dark 
matter in the inner regions of galaxies (see e.g. lTortora et alj2009h 
will make the rejected models even more unlikely, since the analy- 
sis presented here predicts more stellar mass than lensing mass for 
those rejected models. Therefore, our conclusions are robust with 
respect to the complex populations of disk and bulge stars, or to the 
presence of dark matter. 
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